Time-separated entangled light pulses from a single-atom emitter 
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The controlled interaction between a single, trapped, laser-driven atom and the mode of a high- 
finesse optical cavity allows for the generation of temporally separated, entangled light pulses. En- 
tanglement between the photon-number fluctuations of the pulses is created and mediated via the 
atomic center-of-mass motion, which is interfaced with light through the mechanical effect of atom- 
photon interaction. By means of a quantum noise analysis we determine the correlation matrix 
which characterizes the entanglement, as a function of the system parameters. The scheme is feasi- 
ble in experimentally accessible parameter regimes. It may be easily extended to the generation of 
entangled pulses at different frequencies, even at vastly different wavelengths. 

PACS numbers: 42.50.Dv, 32.80.Qk, 32.80.Lg 



I. INTRODUCTION 

The increasing interest in realizing quantum networks 
by means of atoms and photons has risen the issue of 
achieving full coherent control on atom-photon interac- 
tions. Photonic interfaces with atomic ensemble, namely, 
with a macroscopic number of atoms, have been explored 
in several milestone experiments, which demonstrated 
the generation of single photon sources two- mode 
squeezing in the polarization of the emitted light [2|, 
atomic memory for quantum states of light [!, 0,11[ , en- 
tanglement of remote ensembles @, and teleportation 
between light and matter 

Complementary to this approach, photonic interfaces 
using single atoms can take advantage of the large level 
of control that can be achieved on the atomic internal 
and external degrees of freedom. For instance, in mi- 
crowave cavity QED quantum state and entanglement 
engineering using single atoms have been demonstrated 
[1, Q. In the optical regime, milestone experiments 
demonstrated the realization of one-atom laser [l(| [H| , 
revealed the mechanical forces of single photons on sin- 
gle atoms [H, [ID, [l4[ , achieved the controlled interaction 
of a trapped ion and a cavity Il5l llol. yielding single- 
photon generation on demand [171 llq . [l9(, and char- 
acterized the entanglement between a single atom and 
its emitted photon [2(| [2l[ . This latter step was instru- 
mental for establishing entanglement between two distant 
trapped particles by projective measurement of the emit- 
ted photon [22| . Most recently, reversible quantum state 
transfer between light and atoms in a cavity have been 
experimentally demonstrated [23| . These results consti- 
tute relevant progress towards the realization of quantum 
networks with single atoms [24l [25| . 

Possible implementation of quantum networks with 
continuous variables [26j using single atoms as interfaces 
requires the controlled interaction of light with the atom 
external degrees of freedom, which exploits the mechan- 
ical effects of light- atom interactions [13, [H, [29[. Us- 
ing these concepts, in a recent proposal we predicted 



that a single atom, confined inside a resonator, can 
act as a quantum medium, generating quantum corre- 
lations in the emitted light [3(1 [Hj]. Such correlations 
can be two-mode squeezing type of correlations, which 
for bipartite Gaussian systems are synonymous of EPR- 
entanglement [26l . |32| . In particular, under suitable con- 
ditions, two classical light (laser) pulses, temporally sep- 
arated at the input, exhibit two-mode squeezing type 
of correlations at the output of this kind of device, as 
sketched in Fig. [TJ In this case the quantum state of 
the atomic motion serves as intermediate memory which 
mediates the entanglement between the first and second 
pulse at the cavity output. Variation of the laser param- 
eters, driving the atom, allows for tuning the degree of 
entanglement between the pulses. 

In this work, we analyze the efficiency of the pro- 
posal for t emp orally-separated entangled pulses with sin- 
gle atoms |3Ct l3l| by using a quantum Langevin equation 
description. This description permits us to determine 
the correlation matrix and hence the amount of entan- 
glement one obtains using experimentally accessible pa- 
rameter regimes. We show that this proposal is viable to 
existing experiments, hence providing an important step 
towards continuous- variable photonic interfaces with sin- 
gle atoms. We remark that here the atom acts as a source 
of continuous variable "time-bin" entangled pulses, which 
could be an alternative solution for secure communica- 
tion [33j with respect to those employing single-photon 
qubits [3J|. 



This article is organized as follows. In Sec.|TT]the model 
determining the system's dynamics is introduced, and the 
quantum Langevin Equations for the coupled dynamics 
between atomic motion and cavity mode are derived. In 
Sec. IIIII the corresponding correlation matrix for the two 
propagating correlated pulses is determined and the de- 
gree of entanglement characterized. In Sec IIVI entan- 
glement is discussed as a function of the experimental 
parameters. The conclusions and outlooks are presented 
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at the end of the pulse the cavity mode and the atom's 
vibrational motion will be entangled. Cavity decay will 
give rise to a propagating pulse, whose photon number 
is quantum-correlated with the vibrational phonon num- 
ber mm- 

We then assume that at a time T > T\ a second laser 
pulse tuned to the frequency lul — lu c — v drives the 
atom till the time T%- In this regime, one has resonant 
enhancement of the emission of a cavity photon with the 
simultaneous absorption of a vibrational quantum, see 
Fig. LUb). The effective dynamics is described by the 
Hamiltonian 

ff( 2 ) =mx2d t & + H.c, (2) 



FIG. 1: A trapped atom is confined by a harmonic potential 
of frequency v inside a resonator. The atomic dipole is driven 
by two temporally separated laser pulses, whose intensity as a 
function of time is displayed at the bottom, and couple with 
a cavity mode at frequency lu c . The intensities of the two 
emitted pulses at the cavity output are shown as a function 
of time. In this article we show that, in a suitable parameter 
regime, they exhibit quadrature entanglement. 

in Sec. [Vj and the appendices report the details of the 
derivations in Sec. HI] and Sec. IIIII 

II. THEORETICAL DESCRIPTION 

We summarize the basic concepts at the basis of the 
proposal for the creation of pairs of temporally separated, 
continuous- variable (CV) entangled pulses. Let us con- 
sider a single atom in a harmonic trap of frequency v 
and confined inside a resonator, in the setup sketched 
in Fig. [T] Be b, w the annihilation and creation opera- 
tor of an excitation of the quantized motion inside the 
trap. The atomic dipole couples with a cavity mode at 
frequency lo c , which is far-off resonance from the dipole 
frequency uiq- We denote by a, the annihilation and 
creation operator of a cavity photon. A first laser pulse 
at frequency ujl — lo c + v illuminates the atom in the time 
interval [0, T±\. In this regime, the simultaneous emission 
of a cavity photon and of a vibrational quantum of the 
atom is resonantly enhanced, see Fig. [2ja) , and the rele- 
vant dynamics are described the effective Hamiltonian 

frW= iftXioW+H-c, (l) 

where xi is the transition amplitude of the resonant pro- 
cess. This Hamiltonian describes an interaction giving 
rise to two- mode squeezing, i.e., CV entanglement be- 
tween the center-of-mass oscillator and the cavity mode. 
If the pulse is implemented for a sufficiently short time, 
so that photon leakage out of the cavity has negligible 
effects and the dynamics can be assumed to be coherent, 



where \i is the transition amplitude of the resonant pro- 
cess. Again, we assume that cavity decay can be safely 
neglected, so that the dynamics is coherent. When the 
pulse duration 5T = X2 — T is appropriately chosen, the 
quantum state of the center-of-mass motion at time T 
is completely transferred to the cavity mode at time T^. 
Consequently, by cavity decay a second pulse at the cav- 
ity output will be generated, which is entangled with the 

first one [11113 ■ 

This proposal is based on the assumption that during 
the laser pulses the dynamics is essentially described by 
the effective Hamiltonian (fTJ) and ([2]), while detrimental 
effects like atomic spontaneous emission, vacuum opti- 
cal input noise entering the cavity, and fluctuations of 
the trapping potential can be neglected. These assump- 
tions are justified in certain parameter regimes, which 
have been discussed in [3(J EHI • The scope of this work is 
to explore the robustness of the scheme when detrimen- 
tal effects are small but cannot be a priori neglected in 
the dynamical equations. At this purpose, in this sec- 
tion we adopt a quantum Langevin equations treatment, 
taking into account all sources of noise. This permits 
us to determine the correlation matrix for the generated 
pulses, which are discussed in Sec. IIII) and to quantify 
their degree of entanglement for a wide range experimen- 
tal parameters, as shown in Sec HVl 



A. The system 



We consider an atom of mass m, whose center-of-mass 
motion takes place essentially in one-dimension. We as- 
sume, in fact, that the radial potential is sufficiently 
steep, so that the radial motion can be considered frozen 
out. Be the motion along the i-axis, and be x,p the 
position and momentum of the atomic center of mass. 
The center of mass is a harmonic oscillator with angular 
frequency v, whose Hamiltonian reads 

H mBC = fiv(b%+^\ , (3) 
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First pulse Second pulse 

(Ol = tOc + V <D|_ = <D C - V 

FIG. 2: Relevant energy level and resonant transitions dur- 
ing the two pulses. Here, \g) and \e) denote ground and 
excited state of the atomic dipole transition with linewidth 
7, |n c ) the number of cavity photons, and |n m ec) the num- 
ber of vibrational excitations, (a) First pulse: The laser is 
detuned by Ai from the dipole, and by 8\ — v from the 
cavity mode, driving resonantly the two-photon transition 
) — > \g,n mcc + l,n c + 1). These processes gen- 
erate two-mode squeezing between cavity and motion, (b) 
Second pulse: The laser is detuned by A2 = Ai — 2v from 
the dipole, hence driving resonantly the two-photon transition 
\g,n mcc ,n c ) — » \g,n mcc — l,n c + 1). Choosing the pulse du- 
ration properly, quantum-state transfer between motion and 
cavity mode can be achieved, see [27l l2q| . Note that by vary- 
ing the frequency of the second laser pulse and/or the atomic 
transition which it excites, the two output pulses can be at 
different frequencies or wavelengths. 

where x,p are related to the annihilation and creation 
operators b and of a quantum of vibrational energy 
fiv by the relations x — -J%/2mv(b + b') and p = 
i^J Timv / 2{W — b). The relevant atomic internal degrees 
of freedom are the ground state \g) and the excited state 
|e), which form a dipole transition with moment d and 
frequency luq. The transition couples to one optical mode 
of the cavity at frequency oj c and to the laser, a classical 
field, which is a pulsed excitation whose central frequency 
ojl is tuned from a pulse to the next. We define our model 
in the reference frame rotating at wj,, remembering that 
the frequency changes from the first to the second pulse, 
and the two reference frames are hence related by a global 
time-dependent phase. In the reference frame rotating at 
u>l the total Hamiltonian is H = H a + H c + H lnt . Here, 

H a = -hA\e){e\+H mcc (4) 

is the atomic Hamiltonian with A = ojl — <^0j while 

H c = -hSa^a (5) 

describes the dynamics of the cavity mode, with 6 = u>l — 
lu c . The coupling between atom and fields is described 
by 

H iut = h(^B(t)+aB\t)) , (6) 



where a = \g)(e\ and eft = \e)(g\ denote the dipole low- 
ering and raising operators and B(t) is the operator for 
the field degrees of freedom, which we decompose into 
laser and cavity components, B(t) = -Bz(t) + B c , The 
cavity term is B c = g c cos(kxcosO c + </> c )a, where g c is 
the coupling strength and the cavity mode wave vector k 
(k = \k\) forms an angle 8 C with the axis x of the motion. 
The angle <f> c takes into account the position of the trap 
center inside the cavity. From now on we will assume that 
the atomic motion is in the Lamb-Dicke regime, such that 
the atom-photon interactions can be expanded at second 
order in the Lamb-Dicke parameter 77 = k^Jh/2rav. In 
this limit term B c takes the form [3f| 

8 .. fcM .^(i-£. M >«».» + i)) 

—77 cos 6 c g c sin (f) c a(b' +b). (7) 

The laser term is B L (t) = n(t)e ikxcosf>L , where Q(t) is 
the (slowly- varying) Rabi frequency and 9l is the angle 
between the direction of propagation of the laser and the 
trap axis. In the Lamb-Dicke regime this operator reads 

B L (t) = fi(t) ^l-^cos 2 L (26t&+l)j (8) 
+ ir)fl(t) cos 6 L (tf + b). 

B. Quantum Langevin equations 

The full dynamics of the system must take into account 
the coupling with the environment, which is here repre- 
sented by the dipole fluctuations, giving rise to spon- 
taneous emission, by vacuum fluctuations at the cavity 
input, giving rise to cavity decay, and by fluctuations of 
the trapping potentials, which are responsible of damping 
and loss of quantum coherence of the center-of-mass mo- 
tion. Cavity decay is described by the Markovian noise 
operator a m (t), while the heating due to the fluctuations 
of the trap potential is described by the phenomenolog- 
ical Markovian input noise operator b m (t) acting on the 
atomic motion. These two noise sources are mutually 
uncorrelated and have zero mean value, and their only 
nonzero second-order correlation functions are 

(a^ n (t)a in (t'^) = 5(t-t'), (9) 
(b m (t)b m (t')^) = (N + l) 8(t - *'), (10) 
(b m (tyb m (t')) = NS(t-t'), (11) 

where N is mean thermal vibrational number of the ef- 
fective thermal reservoir coupling to the atom center- 
of-mass motion. We denote by Kb the damping rate of 
the vibrational motion. Pure heating corresponds to the 
limit of Kb — > 0, simultaneously with an infinite temper- 
ature of the associated effective reservoir, i.e., N — > 00, 
with 2KbN = Kh, the heating rate, kept constant [37| . 

The input noise terms associated with spontaneous 
emission are more involved because the latter affects both 
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the internal and the motional degree of freedom of the 
atom, due to the presence of recoil. This effect has been 
neglected in the analysis presented in [38[ , which focussed 
on the c.w.-generation of entangled light, and will be 
systematically taken into account in this work. At this 
purpose, the corresponding Langevin force must be de- 
fined. At second order in the Lamb-Dicke parameter 77, 
the Langevin force operator F(t) associated with sponta- 
neous emission is given by (see also its detailed derivation 
in App. A) 



= \H I dcos 



(12) 



X f 1- y cos 2 9{2¥b + 1) + irj cos 9(b^ + b) I , 

where 7 is the spontaneous emission rate of level |e), 
and A/"(cos 9) is the dipole pattern of emission, such that 
J d cos 9Af (cos 9) — 1. The integral is over the angle 
9 between the wave vector of the emitted photon and 
the axis of the motion, and fg n (t) are the zero-mean, 
angle-dependent Langevin forces, with the only nonzero 
correlation function 



(n n (t)fPHt'))=S(9-9')5(t-t'). 



(13) 



It is convenient to isolate the zeroth-order term in the 
Lamb-Dicke parameter in Eq. (fl~2"|) and to rewrite F(t) 
as 



F(t) = r(t)+F nl (t), 

with the zeroth-order term 



f m {t) = / dcos0^(cos0)/j n (i) 



possessing the correlation function 

<f n (t)f nt (*')>= <*(*-*')• 



(14) 



(15) 



(16) 



Term F n i(t) = F(t) — f m (t) is hence at higher order 
in 77. We shall see that spontaneous emission noise is 



essentially described by / m (t), because the higher or- 
der term F n i(t) gives a negligible contribution to the 
Langevin equations in the parameter regime considered 
in this work (see App. B). 

In the frame rotating at the laser frequency the quan- 
tum Heisenberg-Langevin equations (HLE) of the system 
read 1361 



a(t) = (iS - n)a{t) + ia{t) [B\t), a(t)] + V2^a m (t)(l7) 

b(t) = -ivb(t)+ia{t)[B\t)+F\t),b{t)] (18) 
+ia{tf [B(t) + F(t), b(t)} - n b b(t) + V2iTbb m (t), 

a{t)+ia z {t){B{t)+F{t)), (19) 



a(t) 



iA-1 



2ia(t) [B\t) + F\t)] - 2icr t (i) [B(t) + F(t)] 

-7(<r,(t) + l)/2, (20) 
where er 2 = er^rr — 00^ and n is the decay rate of the 
cavity mode. 

In the cases we are going to consider both laser as 
well as cavity mode are tuned far-off resonance from the 
dipole transition, i.e., |A| ^> f2,g, |<5|,7- In this regime 
the atomic internal degrees of freedom can be eliminated 
in second order perturbation theory. Hence, we neglect 
the time evolution of a z , Eq. (|20[) , and approximate 
a z (t) rj —1. Correspondingly, Eq. (fT!))) becomes 

&(t) =-(l~ *A) <t) ~ iB{t) - iF(t), (21) 

whose formal solution is 

tr(t) = e -(*-' ;A )V(0) (22) 

-i / dse-(l- lA ) s [B{t-s) + F(t-s)}. 
Jo 



We insert solution (|22|) into the other HLE and neglect 
the transient term as we are interested in the dynamics 
at times which are much larger than 1/|A|. The resulting 
HLE for cavity and trap oscillators are 



a\t) = -i5a\t)- / ds e~(l +lA ) s [B^ {t - s) + F\t - s)} [B(t),a?(t)] - na^ (t) + V^a" 1 ^ (t) , (23) 
Jo 
pt 

b(t) = -ivb{t)+ / ds e-^- lA >[B{t~s) + F(t-s)} [flt(t)+Ft(t),6(t)] 
Jo 

pt 

- ds e-(l +iA ) s [Bt(f _ s ) + i?t( t _ s )] [ B (t) + F{t),b(t)} - n b b(t) + ^ b b in (t), (24) 
Jo 

I 



where we have not taken care of operator ordering, since, these integral terms will generate only linear contribu- 
as we shall see, within the validity limit of our treatment tions. We now determine the solutions of these equa- 
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tions for the dynamics during the first pulse, between 
the pulses, and during the second pulse. 

1. Dynamics during the first laser pulse 

We first consider the dynamics during the first pulse, 
i.e., in the time interval < t < T%, where we assume 
a square laser pulse with central angular frequency u>li 
and constant Rabi frequency = f^i during this time 
interval. Correspondingly, we denote by A = Ai and 
S = Si the detuning of the laser from the atomic and 
cavity frequency. We assume that the laser is far-off res- 
onance from the atomic transition, i.e., Ai is negative 
and \Ai\ 3> Qi,g,~/. The laser frequency is tuned to the 
value luli ~ uj c + v, namely, the cavity mode is resonant 
with the Stokes motional sideband of the laser light. This 
condition allows to establish a parametric-amplifier type 
of interaction between cavity mode and motion, which 
is selectively enhanced provided that the Stokes side- 
band is spectrally resolved, namely, when vT\ 3> 1 and 
v 3> k, Kb- In particular, we take into account the pres- 
ence of a. c. -Stark shifts Sv, due to the mechanical cou- 
pling with laser and cavity modes, by tuning the laser 



frequency to the value 

wlx = uj c + v' ', (25) 

and hence 5' x = v' , with v' = v + 5v and 8^ = 8\ — 5' , 
and S' accounts for possible a. c. -Stark shifts due to off- 
resonant couplings. The value of 5v is determined in 
a self-consistent way, which is extensively discussed in 
App. [Bl see also |38[. 

Starting from Eqs. (f2"3")l - (|2"4")l . in this parameter regime 
we derive the effective HLE, which describe the coherent 
interaction between the cavity mode and the vibrational 
motion during the first laser pulse, in the presence of 
losses and noise processes due to spontaneous emission, 
cavity decay, and vibrational heating. For later conve- 
nience, we study the equations in the reference frame 
rotating at the cavity-mode frequency, which is obtained 
from the reference frame of the laser frequency by the 
transformation 

a\t) = e i "V(i), (26) 
b(t) = e iv>t b{t). (27) 

In this reference frame, the HLE read 



5 + (i) = X fb(t) - (k + k l ) tf(t) + V2^d m Ht) + V2R L ~a m \t) , 
b(t) = X ia\t) - (K b + K b +1 - K b _j) b(t) + V2ir b b ln (t) + V2K b +1 a™(t) - V^R^o^i) 



(28) 
(29) 



The equations are at second order in the Lamb-Dicke pa- 
rameter, and have been obtained neglecting off-resonant 
terms. The steps of the derivation are reported in App.lBl 
We now define and discuss each term appearing in the 
equations. The effective coupling between motion and 
cavity mode is 

rJig*cos0 c 

Xi = T] (cos#l + itan0 c cos# c ) ,(30) 

Ai 

and corresponds to the Raman processes in which laser 
photons are scattered into the cavity mode with a change 
in the center-of-mass excitation. The new noise operators 
appearing in Eqs. (f2"8")) - (|2"9"j) . a±(t), arc defined as 

a™(<) = f n {t)e- iv '\ (31) 

a l _l l (t) = f m {t)e w '\ (32) 

and describe the coupling with external optical modes 
due to the photons scattered by the atom. They pos- 
sess the same correlation functions of the spontaneous 
emission noise f m (t), and at the time scales of interest, 
v't 3> 1, they are uncorrelated from each other. These 
noise components affect both the cavity mode and the 
vibrational motion. In particular, a™(i) describes quan- 
tum noise associated with incoherent scattering by the 
atom of a cavity photon into the external modes at rate 
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|g c | 2 cos 2 « 



2 7 2 /4+ (Ax 



(33) 



while the corresponding input noise scales with 



g c cos < 



2 7 /2 + i(Ai - v>) ' 



(34) 



where kt, = \kl\ 2 ■ The two noise terms in Eqs. (|31|) - 
(l32|) also affect the atom's motion due to the mechanical 
effects of the scattering of laser photons. The incoherent 
change of vibrational quanta due to scattering of laser 
photons takes place at rates [3!| 



.2 7 



l^ipcos 2 e L 



"±i 



2 7 2 /4+(A 1 ± I /) 2 ' 



(35) 



while the corresponding terms k±i scaling the input noise 
are given by (see App. B) 



R±l = ±Wo 



f^i cos 9 L 



2 7 /2t*(A!±z/) 



(36) 



so that = |k±i| 2 - The formal solution of Eqs. ([28 
at time t = T x reads 
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KTx) 



Ti 



fit (TO - 5 _ 1 (T 1 )at(0)+xI/i(T 1 )6(0)+xt / ds h{T x - s) V& n {s) + V2H b +1 af(s) - V2R\a m \s) 



Ti 



ds g-x{T x - s) [V2K5 int (s) + \/2K L a mt ( s ) 



Xi/i(T x )at(0) +g +1 (T 1 )b(0) + F ds ff+1 (Ti - s) k/2^ m (s) + v^a^s) - V^a™^) 

Jo L 

+Xi / 1 ds h(T x - s) [VW nt (s) + ^2R L ~a m \s) 
Jo 1 



(37) 



(38) 



r 



where we have introduced the time-dependent functions 



-K1S* 



cosh(M) ± ^ sinh(ffit) 



pi 



and the parameters 



«is = 



kid = 



K + K L + K h + K b +1 - K b _ 1 



ID- 



(39) 
(40) 

(41) 
(42) 
(43) 



Equations (|37|) -([38 | describe the dynamics of the cou- 
pled motion and cavity mode during the first pulse in 
presence of quantum noise. By setting all noise terms to 
zero, k, Kj = 0, they reproduce the well-known coherent 
two-mode squeezing dynamics [36| . where entanglement 
monotonically increases as a function of the interaction 
time T\. The presence of quantum noise sets a limit to 
the establishing of these dynamics. 



2. Dynamics between the two laser pulses 

In the time interval T\ < t < T the laser is turned 
off and consequently there are no resonant photon scat- 



tering processes which couple the cavity mode and the 
atom's vibrational motion. The HLE describing the sys- 
tem without laser excitation can be obtained immedi- 
ately from Eqs. by setting = and thus 



Xi 



k_ 1 = 0. The resulting HLE are given by 



h(t) = -(k + k l ) a{t) + V2^a m (t) + \/2R L ht l {t), 

(44) 

b(t) = i5 b b{t) - K b b(t) + V2K~bb in (t), (45) 



where off-resonant coupling between motion and cavity 
mode is neglected, and 5 b is defined in App. |E] This lat- 
ter term is due to the definition of the reference frame 
rotating at frequency u', which compensates the laser- 
induced a. c. -Stark shift on the motion: In absence of 
the laser, this component of the frequency v' is unbal- 
anced. The assumption of neglecting off-resonant cou- 
pling between motion and cavity mode is justified when 
the cavity-mode wavevector is orthogonal to the axis of 
the motion, and thus there is no mechanical coupling. 
In general, it is valid at zero order in the expansion in 
the small parameter rjg cos9 c /\uj c — k>o|, which hence im- 
poses a condition on the integration time T under which 
Eqs. (gl])-([15]) are valid. The solutions of Eqs. ([4l| - ([4"5l) 



a(T) = e- ( - K+KL ^ T - Tl '>a(T 1 )+ / ds e -^ +KL ^ T ~ s ^ V2^a m (s) + V2R* L a"!(s) 
b(T) = e<*i-"»>< r - Tl >&(Ti) + V2^ [ T dse«- K ^ T - s W n (s), 



(46) 
(47) 



r 



where a(Ti) and b{T{) are given by Eqs. (|3"T)) - (|3"8|) . They second laser pulse is switched on. 
yield the observables of interest at time T, before the 
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3. Dynamics during the second laser pulse 

We now consider that in the time interval T < t < T 2 
a square pulse of constant Rabi frequency fl(t) = £1 2 and 
central angular frequency 101,1 illuminates the atom. We 
denote by 62 = UL2 — w Cl an d A2 = ul2 — ^>o the detun- 
ing of the laser frequency from cavity mode and dipole 
transition, respectively. In the limit in which processes 



where absorption of a laser photon and of a phonon is 
resonant with emission of a cavity photon, cj^ 2 ~ w c — v, 
we derive the HLE from Eqs. (fT7 |) -([20 ]) . describing the 
interaction between the cavity mode and the vibrational 
motion during the second laser pulse, in the presence of 
losses and noise processes due to spontaneous emission, 
cavity decay, and vibrational heating, 



a(t) = X 2~b(t) -(k + k l ) 5(t) + V2h~a in (t) + s/2R* L a™(t) , 

b(t) = -x5S(t) - (n b + n b +2 - K b _ 2 ) b(t) + V2JT b b m (t) + V2R b +2 a^(t) - V2R b _ 2 a mj < (t) 



(48) 
(49) 



where parameters X2 and k\ 2 are found from xi an d reported in App. [Bj 

n b ±i by replacing Sli — > Q2, 8\ — > <5 2 , and Ai — > A 2 The solutions of Eqs. (|48]) - (l49p at the end of the second 

in Eqs. ([50)) and (f3"3")l . The details of the derivation are pulse, f — T 2 , read 



a(T a ) - 5 -2(T 2 -r)a(r) + x2/ 2 (r 2 -T)6(r) + x2 

r T 2 



r 2 



/ 2 (T 2 - s) [72^6 m (s) + v^k+^s) - V2R b _ 2 a m \s) 



HT 2 ) 



-xlf2{T 2 - T)5(T) + g + 2(T 2 - T)b(T) + 
rT 2 



-X2 



dsf 2 (T 2 -s) V2K~& in (s) + V2R* L a™(s) 



dsg +2 (T 2 -s) V^bb m (s) + V2R b +2 a^(s) 



V2R b _ 2 ~^ 



(50) 
(51) 



r 



where we have introduced 



9±2(t) = 
fi(t) = 



-K2S* 



-«2S* 



cos(0 2 t) ± ^ sin(0 2 i) 

#2 



1 



sin(# 2 t), 



and 



«2S 



K L + K b 



v +2 



K 2 D 



K + KL - Kb — K +2 



H = \I\X2\ 2 -K 2 2D . 



(52) 
(53) 

(54) 

(55) 
(56) 



Equations (|50| -(|5i p give the cavity mode and motion 
at the end of the second pulse. Setting all decay and 
noise sources to zero, we recover the ideal polariton dy- 
namics during the second pulse, namely a periodic dy- 
namics at frequency \xo\- Ideally, then, the states of 
the motion and of the cavity mode are swapped when 



T 2 -T = ST§, such that |x 2 |<5T 2 ° = (21 + l)n/2, with I 
integer number. At these values the function f 2 reaches 
its maximum, / 2 (£T 2 )\ K .j— = 1, while g_ 2 vanishes, 
3_ 2 (<5T 2 )\ K -.j=o = 0. The effect of decay and noise is 
to damp the oscillators, hence to modify the oscillation 
frequency and the behaviour of functions / 2 , <?- 2 . In 
particular, the maximum value of / 2 (i) is always smaller 
than unity, and the maxima of / 2 do not occur at the 
same instants of time in which <?- 2 (£) vanishes. We opti- 
mize the process by setting T 2 — T such that it fulfills the 
condition g_ 2 (T 2 — T) = 0. Denoting by 5T° pt = T 2 -T 
the time interval fulfilling this condition, it satisfies the 
relation 



5T 2 pt = — arctan — — 

&2 K2D 



(57) 
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III. QUANTIFYING THE ENTANGLEMENT 
BETWEEN THE TWO PULSES 



The two time-separated pulses at the cavity output 
can be considered as two independent modes, even if 
they originate from the same intracavity field at different 
times. In fact, the output field a out (t) is related to the 
intracavity field a(t) by the input-output relation [3t| 



(t) = V2na(t) - a m (t) 



(58) 



and it is characterized by the commutation relation 
[a out (t),a out (t'y] = 6{t - £')■ In order to have a Quan- 
tity directly related to the detected held, we dehne the 
following integrated output field over a generic measure- 
ment time T m [40j], 



rt+Trr, 



t (t,T m ) = 



d£ 



\t>) 



(59) 



is, [a° ut (t, T m ), a° ut f (t',T m )} = whenever \t-t'\ > T m . 
The two output pulses we are interested in are there- 
fore those associated with the operators aj ut (Ti, T m ) and 
Sj"*(T2,T m ). In the present scheme the two pulses are 
temporally separated and therefore it is natural to con- 
sider T m < T — Ti, which automatically warrants the 
independence of the two integrated output modes. 

In order to characterize the entanglement between the 
two pulses, one usually considers the amplitude and 
phase quadratures of the two independent modes, which 
are in this case 



X° ut {T ]: T m ) 



P° ut (T 3l T m ) 



a° ut (T J ,T m )+5,° ut {T J ,T m y 

~a out {T J ,T m )~~a out {T ,T m )^ 
iy/2 



(60) 
(61) 



and construct the correlation matrix 



The field operators a,j ut (t,T m ), form a class of dimen- 
sionless bosonic operators, [a° T ut {t, T m ), d^ ut *(£, T m )\ = 
1, and they commute, i.e., they describe independent 
modes, as soon as they do not temporally overlap, that 



/ tOUt tOUt _|_ £OUt £OUt\ 

■trout _ VSfc s; ~t~ gfc / 
Vhl — ^ 



where we have defined the four-dimensional vector 



(62) 



J 



tout,T 



I 



In Eq. (f6"2"]) the averaging corresponds to taking expecta- 
tion values with respect to the initial state of the system 
and the environment. We now proceed in determining 
its elements. 

We use the definition in Eq. (f59"]) . the input-output 
relation of Eq. (f58|) . and the explicit solution for a(t) 
in the two relevant time intervals, Tj < t < Tj + T m , 
j = 1,2, which is given by Eq. (|46|) . and get the following 
expression of the integrated output fields as a function of 
the intracavity fields at the end of the pulses a(Tj) and 

I 



of the input noises: 



( (T,,T m ) = a(T m )a{Tj) + n in {T 3 ,T m ), 



with the factor 



a(T„ 



and the input noise term 



2k [l 



-(k+k £ )T„ 



(63) 



(64) 



n m (T„T m ) 




Tj+T m 



dt 



l (t) 



K 



2 k 



-(K,+K L )(Tj+T m -t) 



K + K L K + K L 
e -(K+ KL )(Tj+T m -t)-l 



(k + kl) 



(65) 



r 



Using Eq. (|63|) in Eq. (|62|) . we find that the correlation 
matrix y out can be decomposed into the sum of three 
contributions, 



V" 



a(T m ) 2 V + V m + V* 



(66) 



where V m is the contribution due to the input noise term 
n m (Tj, T m ), V mlx is the contribution due to the correla- 
tion between the intracavity fields at the end of the pulses 
a(Tj) and the input noises, and V is the correlation ma- 
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trix for the quadratures of the intracavity fields X(Tj) = 
(2(r J -)+a(r J )t)/>/2, P(T/) = -i{a{T j )-~a{Ttf)/^. Its 
elements have the form 



(67) 



noises a m (t) and a™ (t) one finds that V 1 ™ is proportional 
to the 4x4 identity matrix 



V m = - [1 - a{T m f] S K 



(68) 



where we have defined the four-dimensional vector 
£ T = (X(Ti), P(T{), X(T 2 ), P(T 2 )), and are reported in 
App. [C] Using the correlation functions of the input 



J 



while V mtx has only four nonzero terms, which are all 
identical, V^ ix = V^ ix = V££ ix = V$ ix = V mix , with 



Vraix „ (rp \ 
= Ct(J-m) 



-(k+kz,)(T-Ti) 



'- 2 (T 2 - T) 



IT 

— it 



K + K L 
I 



(69) 



Equation (|69[) shows that the contributions due to corre- 
lations between intracavity fields at the end of the pulse 
and input noise are zero as soon as we choose the optimal 
transfer condition g- 2 {T 2 - T) = of Eq. (JHTJ) (They are 
in any case negligible when (k + kl)(T — Ti) S> 1). 

Let us now consider what is the optimal integration 
time T m , such that V out « V, i.e. the correlation matrix 
at the cavity output reproduces the correlation matrix 
between the intracavity fields. From Eqs. (prJ)) and (|6"5|) 
we see that T m must be chosen such that the quantity 
a(T m ) is as close as possible to unity. In this case, most of 
the intracavity field is detected at the cavity output at the 
end of the pulse, and at the same time the contribution 
of the input noise is negligible. From Eq. (j6~4"|) one gets 
that a(T m ) < a max , where a max ~ 0.9/^1 + k l /k < 1, 
and it is achieved for ~ 1.25/(k + k^). 

In order to establish the conditions under which the 
two output pulses are entangled we consider the loga- 
rithmic negativity E^/, a quantity which has been al- 
ready proposed as a measure of entanglement [4 ill. In 
the continuous variable case Etf can be defined as [43 



where 



E_\f = max[0, — In 2rj 



Y,{V out ) - [T,{V out ) 2 - 4detF ou *] 



1/2 



(70) 



1/2 



(71) 

with Y,(V out ) ee det A+det B-2 det C, and we have used 
the 2x2 block form of the correlation matrix 



-irOUt _ I A C 

1 1 C T B 



(72) 



Therefore, a Gaussian state is entangled if and only if 
r/~ < 1/2. This is equivalent to Simon's necessary and 
sufficient entanglement criterion for Gaussian states of a 
non-positive partial transpose [43j , which can be written 
as AdetV < S - 1/4. 

We finally comment on our choice to quantify the en- 
tanglement in terms of the logarithmic negativity instead 



of EPR variances [32} . The latter would seem a nat- 
ural choice, but in fact they provide an unambiguous 
characterization of entanglement only for simple exam- 
ples of CV two-mode entangled states, such as the two- 
mode squeezed state. Any reasonable entanglement mea- 
sure has to be invariant under local transformations of 
the quadratures of each mode separately. In our sys- 
tem, however, the common definition of EPR correla- 
tions [43, W3| reads 



£ EPR = - [A(X out (T u T m ) - X out (T 2 ,T m )) 2 
+A(P out (T 1 ,T m ) + P out (T 2 ,T m )) 2 ] , 



with A (A) 2 the variance of A, and depends upon the 
chosen set of quadratures, i.e. it does not possess such 
invariance. The light pulses at the cavity output are in a 
two-mode squeezed state only in the ideal limit when the 
noise contributions to the dynamics are negligible, and 
when the ion's motional state created by the first pulse is 
perfectly transferred to the second pulse. Under realistic 
conditions the state of the two pulses is rather different 
from a two-mode squeezed state; therefore using EPR 
variances would give an ambiguous quantification of the 
generated entanglement. Nonetheless, there is a connec- 
tion between squeezing and the measure of entanglement 
provided by the logarithmic negativity. To be more spe- 
cific, a simple and direct quantitative connection between 
EPR variances and E^/ can be found in the case of sym- 
metric bipartite states, i.e., states which are invariant 
under exchange of the two modes (the two pulses in our 
case) [Hf. In such a case, the quantity 2r/~ = exp(—E_\r) 
gives the largest amount of EPR correlations, that is, the 
minimum achievable value of £,epr which can be attained 
in the CV bipartite state by means of local operations, 
i.e. by considering all possible linear combinations of the 
quadratures of each pulse. 



10 



IV. RESULTS 

We now analyze the basic requirements and the ef- 
ficiency of this scheme using parameters accessible in 
present experiments with atoms in resonators. Our con- 
siderations follow and extend the corresponding discus- 
sion in 31]. For comparison, we will discuss along with 
some realistic sets of parameters an idealised case where 
the cavity decay rate is set to a very small value, and the 
noise terms are suppressed. 

The atom's internal degrees of freedom need to pro- 
vide an optical dipole transition which couples to both 
the laser and the cavity mode. A suitable example would 
be an F — «-> F' = 1 closed atomic transition with the 
quantization axis B along the cavity axis, and B, k L , and 
laser polarization El mutually orthogonal. Ideal candi- 
dates would then be alkali-earth-metal atoms or alkali- 
earth-metal-like ions, but other geometries can be found 
which allow for the realization of this scheme using also 
alkali atoms or alkali-like ions. We consider an alkali- 
earth atom in the geometrical configuration 8 L = 0, 
6 C = it/2, and <f> c = 0. This means that the trap cen- 
ter coincides with an antinode of the cavity mode, and 
that the motion takes place along the direction of the 
laser beam and orthogonally to the cavity axis, such that 
there is no mechanical effect of the cavity field on the 
motion. This assumption is not strictly necessary but it 
is made here in order to simplify the discussion. 

In order to favor motional Raman transitions over res- 
onant scattering, we assume a value of the atom-field 
detuning A which exceeds the values of the coupling 
strengths SI and g c \pn, (n is the average number of cavity 
photons), and of the atomic transition linewidth 7. Tak- 
ing a typical value of the linewidth of an optical dipole 
transition, 7 = 2tt x 5 MHz, we choose |A| = 2n x 120 
MHz, Q = 2tt x 10 MHz and g c = 2n x 1 MHz, which 
are accessible values for state-of-the-art experiments with 
trapped atoms or ions in resonators [ljl, [T^, |44[ • 

We also consider the motion to be restricted to the 
Lamb-Dicke regime, with a Lamb-Dicke parameter r\ — 
0.1. This leads to |xi|, \X2\ — 2tt x 8 kHz for the coupling 
constants, while the loss rates associated with the various 
scattering processes are kl ~ k±j — 2n x 0.2 kHz. The 
trap frequency can be set to v ~ 2-n x 1 MHz, which is 
typical in ion trap experiments [ilj . The heating rate of 
the vibrational motion may be estimated as Kh — 2n x 
20 Hz [H, 

The value of the cavity decay rate k must warrant 
the coherent creation of correlations during the pulses 
(nTi t 2 <§; 1), as well as spectral resolution of the side- 
bands (k <C v)- In the examples we discuss we assume 
values that range from k ~ 2n x 20 kHz, which is experi- 
mentally accessible [H, EH, down to 2tt x 1 kHz, which is 
more difficult to reach with present-day technology, but 
serves as an idealised case for comparison. 

Given the parameters, finally the laser pulse durations 
Xi and T2 — T, as well as their separation T must be ad- 
justed in order to (i) create significant entanglement be- 



tween the first pulse and the motion, and (ii) efficiently 
realise the quantum state transfer between the motion 
and the second pulses. At the same time the motion, 
which acts as intermediate quantum memory, needs to 
remain coherent during laser excitation and cavity out- 
put. 

We first focus on the dependence upon the duration 
of the second pulse, T2 — T. The optimization of the 
state transfer from the center-of-mass motion to the cav- 
ity mode discussed in Sec. IIB.3 yielded the condition of 
Eq. {57]), T 2 ~T — ST° pt . Therefore we expect the en- 
tanglement to be maximum around this condition. This 
is confirmed by Fig. [3l where the logarithmic negativity 
is plotted versus the normalized duration of the second 
pulse, (T2 — T)/5T2° pt at various values of the cavity de- 
cay rate. The duration of the first pulse has been fixed at 
T\ = 40 /isec ~ 2/ 1 x 1 1 1 ; while the time interval between 
the two pulses has been chosen to be related to the cavity 
decay time according to T — T\ = 2/(k + kl)', the other 
parameter values are those discussed above. It is evident 
that the entanglement between the pulses is optimized 
when the duration of the second pulse satisfies Eq. (|5Tj) . 




0.5 1 1.5 2 

(T 2 -T)/8T 2 P t 

FIG. 3: Logarithmic negativity Ex versus normalized dura- 
tion of the second pulse (T 2 ~T)/ST° pt (see Eq. ([BT)l) at four 
different values of the cavity decay rate. From bottom to top, 
k — 2tv x (16, 11, 6.4, 0.8, 0.8) kHz; the top solid curve is plot- 
ted for comparison and corresponds to the smallest value of 
ft with all other noise terms set to zero. The other parame- 
ter values are (\A\,j,£l,g c ,u) = 2tt x (120,5,10,1,1) MHz 
and 77 = 0.1, yielding (|xi|, |xa|, «£> «±j) = 2ir x 
(8.5, 8.5, 0.18, 0.18) kHz; the heating rate is n h = 2tt x 20 Hz. 
The duration of the first pulse is fixed at Ti = 40 /js 
~ 2/|xi| _1 , while for each curve, the time interval between 
the two pulses is related to the cavity decay time according 
to T - Ti = 2/(k + k l ) = (20, 28, 49, 327, 400) fxs (bottom to 
top). 



The dependence of the logarithmic negativity £V upon 
the time separation T — T± between the two pulses is 
shown in Fig. 0] Parameter values are the same as in 
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Fig. [3] except that we have fixed the duration of the sec- 
ond pulse at the optimal value T 2 — T = <5T 2 op ' given by 
Eq. (|57|) . Entanglement (i.e. Em) decays essentially lin- 
early, and notably its lifetime is independent from the 
cavity decay rate k. In fact, Em always vanishes when 
T — T\ ~ k^ 1 , i.e. the two pulses are entangled provided 
that their separation is not larger than the vibrational 
heating time . This is not surprising because the 
vibrational motion acts as the continuous- variable quan- 
tum memory mediating the entanglement, and thus the 
heating time limits the achievable coherent storage time. 



1.2 
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[T-TiJOisec) 

FIG. 4: Logarithmic negativity Em versus the time separation 
between the two pulses T — T\ at the same values of the cavity 
decay rate k as in Fig. [3] The other parameter values are also 
the same as in Fig. [3] except that the duration of the second 
pulse has been fixed at the optimal value Ti — T — 6T% P = 
(19, 22, 24, 29, 29) ^s (bottom to top). 



Finally, the dependence of the logarithmic negativity 
upon the duration of the first pulse Ti is shown in Fig. [5] 
The other two timing parameters are, for each curve, 
T 2 -T = 5T° pt as in Fig. H and T - T x = 2/(k + k l ) 
as in Fig. [3] The rest of the parameters have the same 
values as before. Em is always increasing and then tends 
to saturate at a value that, as expected, is larger for 
smaller cavity decay rates. In presence of noise one finds 
an empirical expression for the asymptotic logarithmic 
negativity given by E^ ym ~ In [4|xi|/k] /4. This be- 
havior can be intuitively explained by the fact that the 
first laser pulse entangles the cavity mode and the vi- 
brational motion like in a parametric amplifier, and this 
continuous variable entanglement increases for increasing 
Ti. Cavity losses, however, limit the entanglement gen- 
eration and are ultimately responsible for the saturation 
of the entanglement at larger T\. 



From Fig. [5] it can be seen that choosing a large T x , 
in the saturation regime, has two advantages: i) entan- 
glement is maximized; ii) the scheme is insensitive to 
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FIG. 5: Logarithmic negativity Em versus duration of the 
first laser pulse T\ , for the same values of k as in the previous 
figures. The time interval between the two pulses and the 
duration of the second pulse are set for each curve according 
to T - Ti = 2/(k + Kl) and T 2 - T = ST° pt . The other 
parameter values are the same as before. 

fluctuations of Ti. In must be kept in mind, however, 
that we are considering the entanglement between two 
output light pulses which are counted for a time inter- 
val T^ = 1.25/(/c + kl) starting only when each exciting 
laser pulse has finished (see Sec. III). This means that 
for large values of Ti the detection of the two entangled 
pulses is more difficult, especially when the cavity decay 
rate is large. In these latter cases, the number of pho- 
tons leaving the cavity during the first excitation pulse 
is larger than the average number of photons in the two 
detected output pulses. In the saturation regime T± = 40 
/xsec, the average number of photons per entangled pulse 
is n(Ti) = (at(Ti)a(Ti)) = (12.38,3.65,1.55,0.76) for 
k = 2tt x (0.8, 6.4, 11, 16) kHz, where we have set T* = 
1.25/(k + kl), while the other parameters are as in the 
figures. 



V. CONCLUSIONS 

To conclude, we have characterized the quantum cor- 
relations of two temporally separated entangled light 
pulses, emitted from a single atom inside an optical cav- 
ity, in the set-up first proposed in [30l [3lT | . By means of 
a quantum noise analysis we have quantified the amount 
of entanglement one can extract from this system for ex- 
perimentally accessible parameter regimes, and we have 
shown that the quantum motion of a single trapped par- 
ticle is an efficient quantum medium which creates and 
mediates entanglement on demand between subsequent 
radiation pulses. 

This scheme offers promising perspectives for atom- 
photon interfaces and for devising new cryptographic 
schemes exploiting time-correlated pulses and continu- 
ous alphabets, thus extending those based on time-bin 
entangled photon pairs [HJ]. The scheme may also be 
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easily generalised to the generation of two pulses of dif- 
ferent frequencies, or even at vastly different wavelengths, 
by varying appropriately the frequency or wavelength of 
the second excitation pulse. 

An interesting outlook is to study the scalability of 
the scheme when the number of atoms composing the 
quantum medium is increased in a controlled way, hence 
characterizing which resources the collective excitations 
of the medium may offer for creating entangled light. 



VI. ACKNOWLEDGEMENTS 

The authors acknowledge discussions with Stefano 
Mancini. This work was partly supported by the Eu- 
ropean Commission ("CONQUEST", MRTN-CT-2003- 
505089; "EM ALT', MRTN-CT-2006-035369; "SCALA", 
Contract No. 015714), by the Spanish Ministerio de Ed- 
ucation y Ciencia ("QOIT", Consolider Ingenio 2010 
CSD2006-00019; "LACSMY", FIS2004-05830; "QLIQS", 
FIS2005-08257; " QNLP" , FIS2007-66944; Ramon-y- 
Cajal program), and by the Italian Ministero dell Uni- 
versita e Ricerca (PRIN-2005024254). 



APPENDIX A: MECHANICAL EFFECTS IN THE 
LANGEVIN FORCE 

We derive the Langevin force considering the mechan- 
ical effects of the spontaneously emitted photon on the 
atom. At this purpose, we consider a simple model, con- 
stituted by a dipole of frequency luq at position x and 
Hamiltonian of the center of mass Hcm , 

flat = huj \e)(e\ + H C m 

The dipole is coupled to the modes of the e.m. -field at 
frequency lu s and Hamiltonian 

H em { = ^2 huj s ala s , 

s 

with interaction Hamiltonian 

H int = h^2g s e ik ° x a^a s + H.c. 

s 

where g s is the vacuum Rabi frequency. The formal so- 
lution of the Heisenberg equation for a s gives 

a s (t) = e- iuJst a s (0) - ig s [ dre~ lu)7 'e-' lksX{t ~ T) a{t - t) 
Jo 

(Al) 

We substitute this result into the Heisenberg equation for 
cr, and obtain 



(A2) 



We make now the Markov approximation, assuming that 
the characteristic frequencies of the center of mass mo- 
tion are much smaller than the optical frequencies of the 
e.m. -field, which couple quasi-resonantly with the dipole. 
Hence, in the integral we approximate x(t — r) ps x(t). 
Taking alt - r) w e l "° T a(t), we obtain 

a = -iuj cr + (J z i^2g s e iksX e- luJst a s (0) (A3) 

S 

-a{t)Y,g 2 s e^-^dr 

The integral gives a real term, the linewidth 7 — 
27r^ s gl5{ujQ — u> s ), and an imaginary part, the Cauchy 
principal value shifting the transition frequency. Includ- 
ing this shift in the value of luq wc obtain 



a = -iu a - -a + ia z F (t) 



(A4) 



with the noise source 

F (t)=Y,9se tk ° x e-^ t a s (0) (A5) 

s 

such that (Fo(t)F(,(t r )^) = 7<5(i — £'). Let us now investi- 
gate the form of Fo(t) when the center of mass motion is a 
harmonic oscillator, Hqm = H mec given by Eq. ([3]). Ap- 
plying the Lamb-Dicke expansion, Fo(t) takes the form 

F (t) = Y,9.*- iu ' t a.(P) (A6) 

S 

x (l - y cos 2 6 s (2tfb+ 1) + ir\ cos 9 s (tf + b) I 

with k s x — j] cos 8 s (b' +b), where 9 S is the angle between 
the wave vector of the emitted photon and the axis of 
the motion. Taking the continuous limit of the sum, we 
separate the integrals over the modulus and the polar 
angle, and over the azimuthal angle 9 = 8 S , obtaining 
F (t) fa F(t), with 



F(t) = V7 / dcoB0y/Ar{cosO)fi n (t) 



(A7) 



x I 1- y cos 2 ^(26 t 6 + l) + ir]cos9(b^ + b) ) , 

where Af(cos9) is the dipole pattern of emission, 
/ dcos^A/'(cos6') = 1 and fg n {t) are the angle-dependent 
Langevin forces, 

(ft(t)fPHt')) = 5(6-6')5(t-t'), (A8) 

and which have zero mean value. 



Jo 



^ik s [x{i) — x(t—r)]—iuj s T 



e"*''V(i - t) 



APPENDIX B: EFFECTIVE QUANTUM 
LANGEVIN EQUATIONS 

First laser pulse. Starting from Eqs. ([2"3" ]) -([2"3 ]) . 

which are defined in the reference frame of the laser, we 
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move to a frame rotating at the effective vibrational an- 
gular frequency v 1 ~ v, and we neglect all the terms 
oscillating at v' or larger. This approximation is justi- 
fied in the regime we consider, where we assume that the 
Stokes and anti-Stokes sidebands are spectrally resolved. 



The operators in this reference frame are connected to 
the ones in the reference frame of the laser by the trans- 
formation fit(i) = e iv ' l a^{t), b(t) = e i!/ '*&(t), and their 
equations of motion have the explicit form 



a{t) =i(y' - h) ~a\t) - na\t) + V2k~a in (tf ~ / dse"(*+* Al ) s [B\t - s) + F\t - s)] e w ' 



. cos 4> c (l - 2L cos 2 9 c (2p~b + l))-r)g c sin (f> c cos 6 C (b(t) e - w,t + 6 f {t)e iv>t ^ 

b{t) = i {u 1 - v) b{t) - K b b(t) + V2^b m {t) + / dse-(i-* Al ) s [B(t - s) + F(t - a)] e w>t 

Jo 

irjttl cos9 L +r ) g*sm(f> c cos9 c & ji (t)e- tl/ ' t +T 1 2 b(t)e- W ' t (fll cos 2 9 L + g* cos</) c cos 2 6^5* (i)^ 4 "'* 



(Bl) 
(B2) 



+»/V7(/r t (0 ^s9) e + rfSybWe-^iffHt) cos 2 9) e 

JO 

-W7Q1 cos 9 L + r\g c sin C cos 9 c a(t)e LV ' 1 + ifb(t)e~ w,t (fli cos 2 9 L + g c cos <f> c cos 2 c a(f )e ll/ * ^ 
-»?V7(/rW cos0) e + r, 2 ^jb{t)e- w,t {f^{t) cos 2 0)/ 

I 



d se -U +lAl )s [flt( 4 _ s ) + F \ t - s)] e lu ' 1 



where (. . .)g in the equation for b(t) denotes the aver- 
age over the azimuthal angle 9 with weight given by 
the dipole pattern of emission A/"(cos 9) , and we have in- 
troduced the noise operators a m (t) = e~ lvt a m {t) and 



b m (t) 



Aw tlAn 



b m (t), which are still delta-correlated. We 



use the explicit expression for B(t — s), thereby neglect- 
ing the terms oscillating at v' or faster, and perform the 



time integrals by making the Markovian approximation 
exp{-(7/2 ± iA x + imv')s} w S(s)/(j/2 ± iAi + imv'), 
for m = — 1, 0, 1. 

After long, but straightforward calculations we get the 
final, effective HLE at leading order in the Lamb-Dickc 
parameter, which read 



~a\t) 
b(t) 



1 (z/ - h + 6') tf(t) + x*b(t) - (« + ^l) ~a\t) + V2^~a m (rf + V2R L ~a m {tf + F a , (B3) 
i (1/ - v - Si) b(t) + xia+O) - (n b + K b +1 - k\) b{t) (B4) 
+V2k~bb m {t) + V2R h +1 af(t) - V2K b _x^{t)^ + F b . 



Here, the coefficients 
Xi = 77^13* cos </> 



Xi = V^i9c cos < 



cos 9l 



i 7 /2 



cos 9l 



Ai 



itan< 



• cost 



i 7 /2 



(B5) 



itan< 



. cost 



vy/2 Ai + i 7 /2 



(B6) 



correspond to the Raman processes in which laser pho- 
tons are scattered into the cavity mode with a change in 
the center-of-mass excitation. Since we assume 7 <C |Ai| 



(and also v <C |Ai|), we shall take xi = Xi) an d a P~ 
proximate xi with Eq. ([317)) . The noise operators a±(t) 
are defined in Eqs. (f3"T]) - (|3"2"l) , while the rates Kl, k±! in 
Eqs. (|3"3" ]) -([3"5 ]l . The noise scaling factors R b ±1 of Eq. J3B] 
are generally given by 



Oi cos 9l 



fii(cos< 



v +i 



k_! = -vn X j- 



7/2 - j(Ai + v') 7/2 + iA x 
Oicos^l fi*(cos#' 



j/2 + i(Ai-u') 7/2 -iAi 



but they reduce to the expression of Eq. ([3l>]) because the 
average over the dipole pattern gives (cos9)g = 0. 



14 



The operators F a and Ft in Eqs. (|B3|) - (|B4|) represent 
non-linear noise terms, associated with incoherent scat- 
tering processes. They give rise to a. c. -Stark shift and 
losses, whose effect is in general detrimental for the effec- 
tiveness of the two-mode squeezing processes, and are of 
the form F a ~ rj 2 g 2 a&b/ (A + ij/2), F b ~ rj 2 g 2 al ab / (A + 
17/2), hermitian conjugates, and corresponding input 
noise operators. These terms can be neglected in compar- 
ison with the laser induced Raman scattering processes 
whenever the inequality fix cos 9l g c Vncos9 c is satis- 
fied (n is the average number of cavity photons), which is 
a condition on the parameters and on the geometry of the 
setup. We will assume this regime, and these terms will 
be neglected during the laser pulses. In this regime, the 
frequency shifts of the cavity mode and of the vibrational 
motion read 



6' 



(A 1 - V ')\g c \' 



7 2 /4+(A 1 - i /) 2 
<Sj = 2T7 2 |f2i| 2 cos 2 L Ai 

/ 7 2 /4 + A 2 -!/ 2 1 

l(7 2 /4+A 2 -j/' 2 ) 2 +^' 2 7 2 A ? + 7 2 /4 



(B7) 
(B8) 



These shifts must be taken into account when tuning the 
frequency of the first laser pulse. In particular, since the 
dynamics we seek relies on the resonant two-photon pro- 
cesses, where a laser photon is absorbed and a the cavity 
photon and a vibrational phonon are emitted, hence the 
cavity mode frequency must be exactly at resonance with 
the Stokes sideband of the driving laser, i.e., 



+ 5' + v + 5v. 



(B9) 



Here, Sv can be extracted from Eq. (|B8[) when \6v\ <C 
v, which is satisfied provided that ?7|f2i| <C |Ai| and 
?/|£li/Ai| <C v, and reads 



2A 1 r ? 2 |r2 1 | 2 cos 2 9 L ( 7 2 /4 + A 2 - v 2 ) 



( 7 2 /4- 



Af 



v 2 Y + v 2 ^ 2 



- 27f|f2i| 2 cos 2 6 L 



A! 



A 2 + 7 2 /4' 



(BIO) 



Consequently, 



(Bll) 



determines the effective frequency of the Stokes sideband. 

We notice that the resonant condition (|B9j) . giving 
the relation uj^i — v' = uj c + 5' , gives us a simple re- 
lation between the observables in laboratory frame and 
the "tilded" observables, which are connected by the re- 
lations 



Kb{t) = e- w,t b(t). 



(B12) 
(B13) 



Second laser pulse. The main difference with the 
treatment of the dynamics during the first pulse is that 
the second pulse is set to a different resonance condition, 
e.g. 



UJL2 =UJ C -V , 



(B14) 



where 



and the difference accounts for the a.c- 



Stark shifts induced by the second laser pulse on cavity 
and motion frequency. When condition (|B14[) is fulfilled, 
the cavity mode is resonant with the anti-Stokes motional 
sideband of the laser light, with angular frequency u)l% + 
v" . Spectral resolution of this resonance is warranted 
when vT 2 3> 1 and v 3> k, n b . 

In the reference frame rotating at the frequency of the 
cavity mode, the HLE read 



a(t) = i [v" + S 2 ) o(i) - «»(*) + V2^d m {t) + ds e -(*- 4A2 > [B(t - s) + F(t - s)] e m " 



-g*cos(f> c ( 1 — — cos 2 6 c (2¥b + 1) ] - //<;/. sin <>,. cos 1 



(fe(i)e- w " t + & t (0e^" t ) 

6(f) = i (y" - v) bit) - K b b(t) + V2k~ b b m (t) + I ds e-(i-^) s [B{t - s) + F(t - s)} e'""' 

Jo 

v n* 2 cose L +V9*siri(l) c cose c a\t)e iu " t + r j 2 b(t)e- iu " t (p* 2 cos 2 9 L + g* cos0 c cos 2 e c d\t)e iv " t 



(B15) 



(B16) 



d s e ~(* +iA2 ) s [B\t -s)+F\t- s)] e™" 1 



+ivVl(fe n Ht) ™*0)e + if ^b{t)e~^ ' {f e n \t) cos 2 9) e 
9l +r)g c sm4> c cos6 c a(i)e~ iu " t + rf i b(t)e~ iv " t (fL 2 cos 2 9 L + g c cos C cos 2 9 c a(t)e~ iu " t 
-ivV7(ft(t) ™s9) e + rf^b{t)e- iv,,t {ff{t) cos 2 9) e 



15 



where the noise operators a m (t) and b m (t) are the same before, we finally get the coupled HLE 
as in Eqs. (|Bip ~ (|B2|) . Following the procedure outlined 



o(t) 
b(t) 



= i{v' + 5 2 - 6') o(t) + X 2~b(t) ~ (« + K L ) o(t) + V2^S m (t) + V2R* L d™(t) + F a , 



(B17) 



Z (z/ — 1/ — # 2 ) — — ( K b + K 



+2 



K b _ 2 ) b(t) + V^~b~b in (t) + V2n b +2 ^(t) - \f2R h _ 2 c?-{tf + $U8) 
I 



Let us now define the coefficients appearing in these equa- 
tions. The quantities kl and 5' are given by Eq. (|33f> and 
(|B7|) . respectively. In fact, together with Rl, they do not 
depend upon the properties of the driving laser; more- 
over the two noise operators d l "(t) and a™(i) are given 
by Eqs. (f3"Tj) and The nonlinear terms i 7 ^ and F;, 

are the same as in Eqs. (|B3|) - (|B4|) . and are negligible as 
we take f2 2 cos6*L ^S> g c cos6 c . The coupling constants 
associated with the Raman scattering processes are 



X2 = V^29c cos ( 



X2 = iiQ, 2 g* c cos ( 



cos Ul 



i tan cf> c cos 9 C 



A 2 + v> + i 7 /2 A 2 + i-y/2 



(B19) 



COS t>L 



i tan r cos t 



A 2 + v' - i7/2 A 2 + i7/2 



(B20) 



As we consider the limit 7 <C |A 2 | we shall take \2 = X2 
from now on. The incoherent emission or absorption of 
a vibrational quantum scales with the rates 



^±2 



7 rf cos? 9 L &i 
2 7 2 /4+(A 2 ±*/) 2 ' 



(B21) 



and in the limit |A 2 | ^> 7,^ the rates scaling the input 
noise read 

7 f2 2 C0S#£ 



K_ 2 = ~iV\ 7; 



2 7 /2 + i(A 2 -0 

T fi 2 COS0L 



K +2-^/ 27/2 _ i(A2 + ^ ) - 

Finally, the frequency shift of the vibrational motion 
reads 



6 _ 2A 2 r? 2 |» 2 | 2 cos 2 9 L ( 7 2 /4 + A|- ^ 
2 ~ (7 2 /4 + A 2 -^' 2 ) 2 + ^ 2 7 2 



(B22) 



r/ 2 |fi a 1 2 cos 2 # L 



2A, 



A 2 + 7 2 /4 



For |A 2 | 3> 7, ?7|f2 2 /A 2 | <C f, we find with good approx- 
imation 

b ^ 2A 2 r 1 2 \n 2 \ 2 co S 2 9 L (7 2 /4 + A 2 -^ 2 ) 
2 ~ ( 7 2 /4 + A 2 -^ 2 ) 2 + ^ 2 7 2 



77 2 |fi 2 | 2 cos 2 9 L 



2A 2 



A 2 +7 2 /4' 



determining, together with Eq. (|B7[) . the resonance con- 
dition for the central frequency of the laser pulse, 



UL2 



6' 



(B24) 



For |A 2 | > v, choosing Q 2 = fi l5 then 5| w (5j, 1/' = 1/ 
and the processes leading to absorption of a phonon 
and emission of a cavity photon are resonantly enhanced 
by choosing the frequency of the second laser pulse at 
lol2 = u> c — v'. We consider this regime, as it simpli- 
fies substantially the calculations. When it is not ful- 
filled, one must consider an accumulated phase, which 
gives simply a total phase shift and hence modifies the 
quadratures exhibiting entanglement. 



APPENDIX C: CALCULATION OF THE 
ELEMENTS OF THE CORRELATION MATRIX 



We derive here the elements of the intracavity correla- 
tion matrix V of Eq. (|57|) . using that the cavity modes 
and the vibrational motion are initially in the vacuum 
state, and the fact that input noise is uncorrelated with 
the cavity mode operators at former times. The elements 
read 



Via = 

V n = V 22 = (a*(T 1 )a{T 1 )) + 1/2 
^33 = <a t (T 2 )5(T 2 )) +Rc{(5t(T 2 ) 2 )} + 1/2 
U 44 = (5t(T 2 )a(T 2 )) - Re {(^(T 2 f)} + 1/2 
U 34 = -Im{{a^(T 2 )a^(T 2 ))} 
V 13 = [(o(Ti)a(T 2 )> + (a(Ti)2t(T 2 )) 

+ (5t(T 1 )5(T 2 )) + (at^Jat^a)) + c.c] /4 
V 2i = - [{a(Ti)a(T 2 )) - (a(T x )at(T 2 )) 

- (at(Ti)a(Ta)) + (a(T x )tat(T 2 )) + c.c] /4 
Vu = -i [<5(Ti)5(T 2 )) - (5(Ti)5+(T 2 )) 

+ (tfiTjafr)) - (a t (Ti)a t (T 2 )) - c.c] /4 

^23 = -* [<a(T x )a(T 2 )) + (a(T!)at(T 2 )) 

- (at(7i)a(T 2 )) - (at(T!)at(T 2 )) - c.c] /4 

(B23) Using Eqs. §7$, gHJl and (J5DJI we obtain 



(CI) 
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(at^HTO) = | Xl | 2 ^| /l (T 1 )| 2 + K + 2|4 1 | 2 )^ Tl d S |/ 1 ( S )| 2 ^ (C2) 

<a t (r 2 )5(r 2 )) = | 5 _ 2 (T 2 - T )|2 e -2(K+ K ,)(T-T 1 ) ^t( Tl )a( Tl )) + | X2 | 2 |/ 2 (t 2 _ T) ^ e -2 Kb{ T- Tl ) ^ {Tl) ~ b{Tl) ^ 

+ | X2 | 2 (« h + 2|^ 2 | 2 ) fj~ T |/ 2 ( S )| 2 + ^ | X2 | 2 |/ 2 (T 2 - T)| 2 (l - e -*»V-*») (C3) 

(a(7\)a(T 2 )) = X2 / 2 (T 2 -T)e^-^)( T - Tl ) (aCWTi)) (C4) 
(at(Ti)o+(T 2 )) = = (a(T!)a(T 2 ))* (C5) 
(at(T 2 )5t(T 2 )) = 2 X2 7 2 (T 2 - T).g_ 2 (T 2 - r) c -K+^+^+«)(^i) ^(^(Tr)) (C6) 
(a(T!)at(T 2 )) = ,g_ 2 (T 2 - T ) e -^+^ T - T ^ (C7) 
<o+(Ti)o(T 2 )> - (a(Ti)5t(T 2 ))* (C8) 

I 

where we have used Kh = Kb(2N + 1), and that termine the above relations as a function of the initial 
Kh 3> Kb- ln deriving these relation we also used that conditions, we need 
<5t(Ti)6(Ti)) = (5(Ti)6t(Ti)) = 0. In order to fully de- 



^(Ti)^)) = \ X i\ 2 \h(Ti)\ 2 + n h j\s |. 9+1 ( S )| 2 + 2«| Xl | 2 ^ Tl d S |A( S )| 2 (C9) 

+2 / ds -Xi«l/i(«)| 
Jo 

o(Ti)6(Ti)) - ,g_ 1 (T 1 ) Xl / 1 (T 1 ) + Xi (k /i - 2n b ) J ' ds h{s)g +1 {s) (CIO) 



+2 ( 1 As {R b * lX ifi(s) - Klg-^s)) (« b _ l5+1 (,s)-S iXl / 1 ( S ))+2« Xl / Vi(*)/i(«) 
Jo Jo 



I 

In deriving these expressions we used that the HLE for long times, when the perturbative expansion loses 
conserve the commutation relations. This is not fulfilled validity. 
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